Time-reversal anomaly and Josephson effect in time-reversal invariant topological 

superconductors 
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Topological superconductors are gapped superconductors with protected Majorana surface/edge 
states on the boundary. In this paper, we study the Josephson coupling between time-reversal 
invariant topological superconductors and s-wave superconductors. The Majorana edge/surface 
states of time-reversal invariant topological superconductors in all physical dimensions 1, 2, 3 have 
a generic topological property which we named as time-reversal anomaly. Due to the time-reversal 
anomaly, the Josephson coupling prefers a nonzero phase difference between topological and trivial 
superconductors. The nontrivial Josesphon coupling leads to a current-flux relation with a half 
period in a SQUID geometry, and also a half period Fraunhofer effect in dimension higher than 
one. We also show that an in-plane magnetic field restores the ordinary Josephson coupling, as a 
sharp signature that the proposed effect is a consequence of the unique time-reversal property of the 
topological edge/surface states. Our proposal provides a general approach to experimentally verify 
whether a superconductor is topological or not. 



In recent years, topological states of matter such as 
topological insulators (TI) and topological superconduc- 
tors (TSC) have attracted tremendous theoretical and 
experimental interest [H-Q ■ The first example of TSC is 
the time-reversal breaking (p + ip)-wave superconductor 
of spinless fermions in two-dimensions More recently, 
a new class of TSC was proposed in three-dimensional 
(3d) time-reversal invariant (TRI) superconductors 0-0], 
which has time-reversal invariant pairing order parame- 
ter in the bulk, and two-dimensional massless Majorana 
fermions with linear dispersion on the boundary. Can- 
didate materials of TSC include the 3 He B phase HrHI, 
Cu-doped Bi 2 Se 3 [l-[ll| and p-type TlBiTe 2 [ll- 

A robust topological state of matter should be char- 
acterized by topological effects that qualitatively distin- 
guishes it from trivial states. The first topological effect 
for TRI TSC was proposed in Ref. Q which considered 
the time-reversal property of a pair of Majorana zero 
modes in a TRI vortex defect of 2D TSC. As a conse- 
quence of the Majorana zero modes, the time-reversal 
symmetry maps a state with odd number of fermions in 
the zero modes to one with even number of fermions. 
This is due to the following equation 



T 1 i7ot7oi7"= -«7ot7cu. 



(1) 



where 7o-|-,4. are the Majorana zero modes which form 
a Kramers pair: T -1 7otT = 704., T -1 704. T = -701- 
1701-704, is the fermion number parity operator with eigen- 
value ±1. The fact that time- reversal transformation 
changes fermion number parity is generic for any Kramers 
pair of Majorana zero modes, which applies to the 
edge/surface states of TRI TSC in all physical dimen- 
sions D = 1,2,3. Since fermion number should be pre- 
served microscopically by time-reversal transformation, 
the anti-commutation between time-reversal and fermion 
number parity given by Eq. ([1]) is an emergent property 
and cannot be modified by any local perturbation to the 



TSC. This property requires the gaplessness of the sur- 
face states, since the lowest energy state with even and 
odd number of fermions are required to be degenerate. 
Conversely, opening a gap at the surface state necessar- 
ily breaks TR symmetry. We name this property of TSC 
as time-reversal anomaly due to its similarity to the par- 
ity anomaly [l3l [T3|. The time-reversal anomaly can in 
principle be measured if one can measures the fermion 
number parity near an edge. However, it is practically 
difficult to measure fermion number parity. Other topo- 
logical effects have been proposed [15l - ll7| but their exper- 
imental observation are also highly nontrivial. 

In this paper, we translate the time-reversal anomaly 
into a topological Josephson effect which is easy to mea- 
sure in current experimental techniques, and can be used 
to distinguish TSC from ordinary SC in all physical spa- 



tial dimensions D = 1,2,3. Similar to Rcfs. [15l |17|. 
we consider the Josephson coupling between TSC and s- 
wave superconductors. The regular Josephson coupling 
in such a junction is very weak due to different pair- 
ing symmetry, but the Majorana surface states induce a 
nontrivial Josephson coupling which is qualitatively dif- 
ferent from the ordinary case. Since the surface state is 
protected by time-reversal symmetry, it remains gapless 
if the relative phase between TSC and s-wave supercon- 
ductor is or 7r, while it generically becomes gapped for 
8 =/= 0,7r. Since opening a gap in the surface states can 
save energy, this simple reasoning leads to the dramatic 
conclusion that the lowest energy state of such Joseph- 
son junction does not occur at relative phase 8 — or n, 
but somewhere between and tt. Furthermore, the time- 
reversal anomaly Eq. (JT]) tells us that the ground state at 
relative phase 8 and —8 have opposite fermion number 
parity, since they are related by time-reversal. Conse- 
quently the system cannot stay at the lowest energy state 
if the fermion number at the junction is conserved. We 
studied the behavior of the topological Josephson cou- 
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pling for both systems with a conserved fermion num- 
ber and those with a fluctuating fermion number which 
may be induced by disorder states in a real system. For 
fixed fermion number the Josephson effect has a phase 
shift determined by the fermion number. For fluctuat- 
ing fermion number the Josephson effect has a doubled 
frequency. In both cases the junction is qualitatively dif- 
ferent from an ordinary Josephson junction, which can 
be detected in the behavior of a SQUID (superconduct- 
ing quantum interference device) and in the Fraunhofer 
effect (for D = 2,3). In particular, in D = 1 or D = 2 
we show that an in-plane magnetic field comparable with 
the pairing gap at the Josephson junction can drive the 
junction back to a normal one, which provides a clear 
signature that time-reversal symmetry plays an essential 
role in the topological Josephson effect. [18] 
Model Hamiltonians 

Although our proposal applies to a generic TSC, for 
the concreteness of our discussion we would like to start 
by defining some model systems for TRI TSC. The TSC 
in D = 1,2,3 can all be realized in the following lattice 
Bogoliubov-de Genne (BdG) Hamiltonian 

D 

Hp =-/i c ™ Crs ~ 1 Y X)( c H-S«s Crs + h - c -) 

rs i— 1 rs 

D 

+ Y}^-^ ■ ^ss'cl+eAs' + h - C "]> ( 2 ) 
2—1 rss' 

In D = 2 (3), the model is defined on a square (cubic) lat- 
tice, so r denotes the lattice sites and e,'s (i = 1, 2, .., D) 
are the orthogonal unit vectors, cr = (cr x , a y , a z ) denotes 
the Pauli matrices, and s =f, J, the electron spin index. 
The TRI TSC phase requires the chemical potential is in 
the range -2Dt < fi< -2{D - 2)t. 
Josephson effect in ID TSC 

We start from ID TSC, which gives us the clearest ex- 
ample of how the Josephson coupling between the TRI 
TSC and the s-wave SC arises from time-reversal sym- 
metry breaking. ID is distinct from 2D and 3D in the 
sense that there is a finite gap separating the Majorana 
zero modes with other quasi-particle excitations. The 
candidate material for the ID TRI TSC include a class 
of quasi-lD organic superconductors [l9| and LiMO [20| . 
As will be discussed below, such a gap makes it possible 
to probe the "TR anomaly" of Eq.([T]), i-e., the fact that 
time-reversal changes the fermion number parity of the 
Majorana zero modes. 

The ID Josephson junction we consider here is de- 
scribed by the Hamiltonian H = H p + H s + Hi, with 
H p the Hamiltonian of a ID TRI TSC chain given by 
Eq ll4l H s that of an s-wave superconducting chain and 
Hi a weak spin-conserving hopping between the first site 
of TSC chain and the last site of s-wave chain: 

Hi = -Y,m) C {j ma +h. c .\. (3) 



For small coupling St we can study the Josephson cou- 
pling by perturbation theory. Since the TSC and the 
s-wave SC are both gapped, the main effect of a small St 
is to couple the Majorana zero modes at the boundary of 
TSC. We start from a special case fj, = 0, A = t for the 
p-wave chain, in which case the TSC Hamiltonian can 
be written as H p = -iE„a( c "+i>" + c i+i,J(v- c L) 
|21[, in which case the Majorana zero modes are com- 
pletely localized at the boundary site: 7o CT = c\ a + c\ a . 
Therefore the interchain hopping Hi can be written as 
Hi = - EAt {70. - (ci.a - c\j}f ma + h.c.]. Only the 
zero mode terms contributes to the second order pertur- 
bation, giving the simple effective Hamiltonian 

H eff = J (*7ot7o|) sin <p. (4) 
i J = (St) 2 J dt (Tf mf (t)f mi (0)) A , =ilAll , 

where A' is the pairing gap of the s-wave chain and (j) its 
phase; our gauge is set to give us a real St and (T . . .) de- 
notes the time-ordered expectation value. This effective 
Hamiltonian has time-reversal symmetry, both sin <fi and 
*7ot7oi being time-reversal odd. 

The effective Hamiltonian of Eq.Q is valid for generic 
parameters fx, A, for which the Majorana zero mode 70. 
is not completely localized at the first site. This is be- 
cause the zero mode term of the inter-chain coupling 
Hi = -uoYJ2<r("/o<jf m a + h.c.) + ... remain identical 
up to a constant factor that comes from the mode ex- 
pansion Cla = "070<t + J2n^o( U nlna + Klta)l where 

Inn {ilia ) ' s are annihilation (creation) operators of finite 
energy quasi-particles. 

Due to its dependence on the local fermion number 
parity «7ot7oj,j the effective Hamiltonian Eq.(j4|) can give 
us a nonzero Josephson current at 4> = 0,7T. This is 
because there is a limit where 170-1-704, is quasi-conserved 
so our Josephson coupling should be taken to be E± = 
±Jsin</> for 170^704. = ±1 respectively. This gives us 
the Josephson current of I((j>) = ±J C cos <f>, where I c = 
2eJ/h, which is nonzero at (j) = 0, tt. This is the physical 
signature of our TR anomaly - the local fermion number 
parity eigenstate breaking time-reversal symmetry. We 
emphasize that this is a generic property of the TRI TSC. 

The above Josephson current from TR anomaly is 
physically observable for a sufficiently large Josephson 
frequency uij — 2eV/Ti. In a clean system, the local 
fermion number parity is absolutely conserved at T = 0, 
while in a realistic system, it may change, e.g. by a 
fermion leaking out of the zero modes to some impurity 
sites. We denote the rate of fermion parity changing pro- 
cess as 1/t. Therefore in the fast limit wj ^> 1/r, 170^704. 
would not change in many Josephson period, i.e. it is 
quasi-conserved. On the other hand, in the slow limit 
ujj <C 1/t, the system would stay at the ground state 
energy Eq = — J|sin0| of the the effective Hamiltonian 
Eq.((3]), as it has enough time to relax into the lowest 
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FIG. 2. The Zeeman field effect on the Josephson current 
in the slow limit for St /A = 0.2. A.: The Zeeman field is 
applied along the junction with the magnitude of gush/ A = 
(black dashed), 0.005 (blue), 0.01 (red), and 0.02 (brown). 
B.: The Zeeman field of magnitude gfiBh/ A = 0.02 applied 
along the junction (red) and perpendicular to the junction 
(blue). While the field along the junction changes the location 
of discontinuity and eventually removes it altogether, the field 
perpendicular to junction has negligible effects. 



FIG. 1. A.: Schematic picture of a DC voltage-biased SQUID. 
B.: Dependence of the current oscillation amplitude on the 
flux C. & D.: Numerical calculation of the DC SQUID 
I/I c as a function of <f> for the flux of $/$o = 0, ^, j, |, \ in 
blue, purple, brown, green, and black respectively; C. is for 
the slow limit and D. for the fast limit. 



energy state. This gives us the current phase relation 



of = J c ^ (— |sin0|) = -7 c ^^cos0, which has 
doubled frequency 2toj with discontinuities at (f> = 0,7r. 

To relate this behavior of Josephson junction to physi- 
cal observable effects, we now consider a DC SQUID built 
from two of such Josephson junctions. We now consider 
the TRI TSC chain to be finite and is coupled to the 
s-wave chain through hopping at the both ends of the 
chain, with flux $ through the SQUID, as shown in Fig. 1 
A. Each of two junctions has a Kramers doublet of Majo- 
rana zero modes 700- and j' Qa , where j' Qa = —i(cN a —c' Na .) 
for the case of /1 = 0, A = t. Therefore, the effective 
Hamiltonian of this SQUID is 

H ef f = -J[N F iSin((t>— 7T$/$ ) +-/V F 2Sim> + 7T$/$o)], 

(5) 

where $0 = h/2e is the unit quantum flux and N F i = 
*7ot7oj,j Np2 = — *7o-[-7o4. arc two local fermion number 
parities. The behaviors of this SQUID in the fast and 
slow limits are qualitatively different. In the slow limit 
ujjT <C 1, the system stays at the lowest energy state with 
the energy Ej = - J[\ sin(</)-7r$/$ )| + | sin(^+7r$/$ )|], 
so that the current-flux relation is 



- = -E 



± I sin 



■ cos 



As shown on the left plot of Fig. I B., there are sharp 
jumps in the current four times every period due to the 
fluctuation of the electron parity of each junction except 
at $ = n<f>o and $ = (n + l/2)$o where the electron 



parities of two junctions fluctuate together to give us the 
maximum current oscillation amplitude. 

Due to the TR anomaly, whether this DC SQUID be- 
haves like a normal SQUID or a 7T-SQUID in the fast 
limit ojjt ^> 1 is determined hysteretically. As we have 
Josephson oscillation without the fermion number par- 
ity fluctuation in the fast limit, we have four possible 
current-phase relations depending on the eigenvalues of 

JV.F1, N F2l 



'F2> 

-2A r F i/ c cos(7T$/$ o )cos0, N F1 N F2 = +l 
-2A^ F i/ c sin(7r$/$ o )sin0, N F iN F2 = -I 



(7) 



We see that we have a normal SQUID for the even total 
fermion number parity N F iN F 2 = +1 and a 7T-SQUID 
with zero supercurrent at integer flux but maximal su- 
percurrent at half-integer flux for the odd total fermion 
number parity N F iN F 2 = — 1. We now note that even in 
the slow limit, the SQUID is in the N F iN F 2 = +1 eigen- 
state at $ = ?i<f>o and the N F iN F 2 = H — I eigenstate at 
$ = (n + l/2)$o as N F i and N F 2 fluctuate together at 
these flux values. Therefore we can obtain the 7T-SQUID 
is by ramping up the DC voltage from the slow to the 
fast limit while the flux is fixed at $ = $o/2, the normal 
SQUID by going through the same process at zero flux. 

The effect of Zeeman field also reveals the unconven- 
tional nature of this Josephson junction. The Zeeman 
field breaks time reversal invariance and provides an 
alternative channel for gapping out the Majorana zero 
modes. To the leading order, our effective Hamiltonian 
of a single junction is modified to 



H eff = (i-yofYoi)(J sine)) + .g^sh • n), 



(8) 



where g is the g-factor of the normal state of TRI TSC, 
[iB the Bohr magneton, and h is the effective Zeeman 
field, and n is a unit vector along the junction 0, Eil ]. 
We can see that there is no level crossing between the 
even and odd parity states once h • n > J/(g/j,B), giv- 
ing us a normal Josephson junction with no Majorana 
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FIG. 3. We have set St/ A = 0.7. A.: Numerical results for the 
current-phase relation I((f>). B.: Numerical results for I(<fi) 
with the Zeeman field g^, B h/A = (black), 0.02 (red), 0.06 
(blue) along the junction. C: The scheme for the Fraunhofer 
diffraction experiment, hout gives us the flux $ through the 
junction, hi n gives rise to Ea. ()12p . D.: The Fraunhofer diffrac- 
tion pattern with the Zeeman field gfish/A = (black), 0.02 
(red), 0.06 (blue) along the junction; the filled curves calcu- 
lated from the current-phase relation Eq. (|13|l and the dots 
to the numerical calculation. 



zero modes. The numerical calculation shown in Fig. 
[2] confirms both the sensitivity to the Zeeman field di- 
rection and the absence of level crossing at sufficiently 
strong field. Lastly, we note that the effect the Zeeman 
field adds to the Josephson coupling imp lies both the 
spin accumulation at the iunction[23l. |24| and the Zee- 
man field induced s-wave pairing on the boundary of the 
TRI TSC. We will show that the same effects exist when 
the Majorana boundary states gets gapped in the higher 
dimensions. 

Generalization to higher dimensions 

In higher dimensions, we can still detect through the 
Josephson effect a Zeeman field induced s-wave pairing at 
the boundary of a TRI TSC. This effect does not come 
from the Majorana zero modes; in fact, in the higher 
dimensions we are always in the slow limit lojt <C 1 as the 
boundary state of of a TRI TSC has macroscopic number 
of modes forming a continuous spectrum. However, the 
Josephson effect can be studied by considering a weak 
tunneling between TSC and an s-wave superconductor 
just as in ID. For 2D the tunneling Hamiltonian can be 

written as Hi = - J2n=i T,A( St ) c U^ n)f a {m, n) + h.c] 
(n = 1,2,..., N y labels the boundary sites), and it has a 
similar form for 3D. In both 2D and 3D, this coupling 
can only open a gap when the relative phase 6 ^ 0, 7t 
just as in ID. 

The period doubling of the Josephson coupling in the 
higher dimension can be obtained for the weak tunneling 
limit by the second order perturbation theory exactly 



analogous to that of the ID case. Here we will omit the 
detail of the perturbation expansion and directly write 
down the effective edge theory induced by the coupling. 
For the ID Majorana edge state of 2D TSC, the effective 
Hamiltonian is 

H c s = ^2\vk (7_ fet 7fc t -7_fe4,7^)+«(Msin0)7_ fe 4.7 fet ] . 
k 

(9) 

As required by TRS, the mass term is proportional to 
sin 6; in fact it is roughly the same as the ID Josephson 
coupling and therefore also proportional to (St) 2 . The 
2D Josephson energy is determined by the ground state 
energy of the effective model ((9]) defined with respect to 
that of the M = system: 



Ej(<t>) 



-^2 -^Jv 2 k 2 + M 2 



sm 



(10) 



In the thermodynamic limit the sum becomes an integral 
and we obtain the energy and current phase relation 



Ej(4>) = 



2tt 



1 + log 



4IAI 



log 



M 2 sin 2 
4IAI 2 



sm 



M 2 sin 2 



sin 20, 



(11) 



where I c cx M 2 . While the current-phase relation does 
not have discontinuity as in ID, it does have a logarith- 
mic singularity at 6 = 0, 7r; this also shows up in the 
numerical calculation plotted in Fig. [3] A. As a conse- 
quence of the contribution from macroscopic number of 
modes, the Josephson current has a doubled frequency. 
Eq. (fT0| remains valid for junctions between 3D TSC 
and s-wave superconductors, but there it gives us a 2D 
integral. Consequently, in the 3D case the junction has a 
doubled frequency without singularity: 1(6) = I c sin 26. 

The in-plane Zeeman field in 2D can erase out this 
period doubling, thus affecting the TSC Josephson junc- 
tion more dramatically than in ID. Analogous to EqJHl 
the Zeeman effect leads to an additional contribution to 
the gap of the edge state in the effective Hamiltonian © 



Hz = ^2 *7-fci7fet9A i s h • 



(12) 



with n the unit vector along the edge normal. This mod- 
ification will induce an ordinary 2-7r periodic Josephson 
coupling, since the current-phase relation is modified to 



log' 



4|A| 2 /M 2 



l( p + h/M) 2 



2h 

sin 26 H cos < 

M 



(13) 

where h = j^h • h. In fact this formula tells us that for 
\h\ ^> \M\ our Josephson coupling becomes essentially 
conventional, with the effective critical current nearly 
proportional to the in-plane field. Given \M\ oc (St) 2 , 
this crossover occurs in the \h\ St regime, which is 
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confirmed by the numerical calculation shown in Fig. [3] 
B. The underlying reason is that the Zeeman energy term 
of Ea. (Tp2")) induces s-wave pairing correlation, which can 
be seen clearly in the real space basis. In fact, for the TRI 
TSC of Eq.® we cannot gap out the Majorana boundary 
state without inducing the s-wave pairing correlation. 

This in-plane Zeeman field induced change in period- 
icity can be also seen in the Fraunhofer diffraction in 
a single Joscphson junction. When the magnetic field 
is completely perpendicular, the Fraunhofer diffraction 
pattern for 2D has zeros at both n$o and (n + l/2)$o- 
We can see this by noting that the 2D Josephson current 
has only even harmonics and therefore can be expanded 
as = X) m / 2m[sin(27rm$/$o)/27rm$/<I>o] sin 2710. 

However if one measures the Fraunhofer pattern in a 
canted magnetic field with both perpendicular and in- 
plane components, Eq. (|13j) gives us nonzero odd harmon- 
ics in the current phase relation as shown in Fig. [3]B., 
which leads to non-zero critical currents at half-integer 
fluxes as shown in Fig|3](d). This restoration of the ordi- 
nary Josephson coupling provides a sharp experimental 
signature of the TSC. 
Summary 

We studied the topological Josephson coupling be- 
tween the TRI TSC and the conventional superconductor 
for D = 1,2,3 which is qualitatively different from or- 
dinary Joscphson effect. The effects proposed originate 
from the time-reversal anomaly-the anticommutation be- 
tween the time reversal operation and the local fermion 
number parity of the topological edge/surface states. In 
ID, the DC SQUID made of two such junctions can hys- 
teretically behave either a normal or a it SQUID. In 2D 
and 3D the Josephson coupling has a ir periodicity, lead- 
ing to a half period Fraunhofer pattern and a half period 
DC SQUID. Since the topological coupling is determined 
by the time-reversal protected edge/surface states, a Zee- 
man magnetic field breaking time-reversal symmetry is 
predicted to restore the normal behavior in the junction, 
and induce a cross-over of the topological Josephson junc- 
tion back to a normal one. 
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Models of TSC in D = 1, 2. 

Here we present the explicit form of the Hamiltonian 
@ for D = 1 and D = 2. For our numerical calculation 
for D = 1 in Figs. 1 and 2, the BdG Hamiltonians for 
our two superconducting chains are 

H p =- fiJ24 s c rs - tJ2(cl +iis c rs + h.c.) 

rs rs 

+ E( Ac r+- c -+h.c), (14) 
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for the TRI TSC and 



H s=- M' E fish's * J2(ti+*J*s + h ' c -) 

r's r's 
+ E( A '/rVr'i+ h - C -)- (15) 



for the s-wave chain; note that Eq. tfH)) can be obtained 
from the D = 1 case of Eq. ([2j) by applying — ir/2 spin 
rotation around the z axis. We have set \t\ = \t'\ = |A| = 
A'|, /i s = fi p — 0, \St/t\ = 0.2, and both chains consisted 
of 100 lattice sites. For the numerical calculation in Fig. 
2, we add the Zeeman field term 



Hz = -g^B E( h • °") 



(16) 

For the D = 2 numerical calculation in Fig. |3l we used 
for the TRI TSC the BdG Hamiltonian 

H P = - 4 s c rs - t E E( c r+ss c ™ + h ' c ') 

rs e=x,y rs 

+ AEK4+x S 4 - ^ r+ys 4) + h-c], (17) 



which can be obtained, just as in D = 1, by applying 
— 7r/2 spin rotation around the z axis. This BdG Hamil- 
tonian gives us p x + ip y pairing between spin up electrons 
and p x — ip y pairing between spin down electrons. The s- 
wave superconductor on the square lattice has the Hamil- 
tonian 

H s =-/i'E fr'Jr's ~ t' E E (/r'+es/r's + M 
r's r's e— x.y 

+ E(A7 r %/ r V+h.c). (18) 



The parameters we used were |A| = |A'|, |t/A| = 
|i'/A'| = 3, /i/|A| = m'/|A'I = -2, and |<yt/t| = 0.7. 
We find that these parameters open up the edge state 
gap of il//|A| = 0.016 at = ?r/2; the fact that this de- 
viates within an order of magnitude from the crossover 



value of h indicated by Fig. [3]B. is because Eq. (fT"3|) as- 
sumes that of the edge state is linear and that the mass 
term is fc-independent, neither of which is strictly true. 
For obtaining the result shown in Fig. [3] A. and B., we 
have set our superconductors on the lattices of 100 sites 
along the junction (the y-direction) and 30 sites perpen- 
dicular to the junction (the ^-direction) , and imposed a 
periodic boundary condition in the y-direction. For the 
numerical calculation in Fig. [3] D., we have 10 lattice 
sites along the y-direction (without periodic boundary 
condition) and 12 sites along the x-direction; in addition 
the perpendicular flux $ modifies the hopping between 
the two superconductors: 



H x = -EE^) 6 IT ° cl(l,n)U(m,n) + h.c.]. 



n—1 a 



(19) 

For the analytic calculation in Fig. [3]D., we used Eq. 
([T3| with M = 0.02/sin(47r/25), which is the value we 
infer from the numerical results for Fig. [3]B. 



Additional contributions from ordinary Josephson 
effect. 



It should be noticed that ordinary, bulk-to-bulk 
Josephson coupling between TSC and s-wave SC is 
possible through the spin-orbit coupled hopping across 
the junction, which would take the form Hsoc = 
EmxE s =±i \[crcl(l,n)f s {m,n + s) + h.c] in D = 2; 
this would lead to an ordinary contribution I(4>) oc sin <fi. 
From perturbation theory, this conventional Josephson 
coupling will be much smaller than the second order 
Josephson coupling through the Majorana edge state in 
the limit of \X/5t\ <C |(5£/A| 2 ; this is consistent with the 
singlet-triplet Josephson junction studied by Asano et 
al. [26( We have numerically studied the interplay between 
topological and regular Josephson effects with the spin- 
orbit coupling strength set at A/|A| = 0.2, which surely 
overestimates A, with other parameters set the same as 
it was for Fig 3 (a) and we found the first harmonics to 
be less than 20% of the second harmonics, even though 
\X/8t\ and |<5t/A| 2 are of the same order of magnitude 
with our parameters. 



